Refining the evolutionary tree of the horse Y chromosome.

The Y chromosome carries information about the demography of paternal lineages, and thus, can prove invaluable for retracing both the evolutionary trajectory of wild animals and the breeding history of domesticates. In horses, the Y chromosome shows a limited, but highly informative, sequence diversity, supporting the increasing breeding influence of Oriental lineages during the last 1500 years. Here, we augment the primary horse Y-phylogeny, which is currently mainly based on modern horse breeds of economic interest, with haplotypes (HT) segregating in remote horse populations around the world. We analyze target enriched sequencing data of 5 Mb of the Y chromosome from 76 domestic males, together with 89 whole genome sequenced domestic males and five Przewalski's horses from previous studies. The resulting phylogeny comprises 153 HTs defined by 2966 variants and offers unprecedented resolution into the history of horse paternal lineages. It reveals the presence of a remarkable number of previously unknown haplogroups in Mongolian horses and insular populations. Phylogenetic placement of HTs retrieved from 163 archaeological specimens further indicates that most of the present-day Y-chromosomal variation evolved after the domestication process that started around 4200 years ago in the Western Eurasian steppes. Our comprehensive phylogeny significantly reduces ascertainment bias and constitutes a robust evolutionary framework for analyzing horse population dynamics and diversity.

Y Chromosome Haplotypes Enlighten Origin, Influence, and Breeding History of North African Barb Horses.

In horses, demographic patterns are complex due to historical migrations and eventful breeding histories. Particularly puzzling is the ancestry of the North African horse, a founding horse breed, shaped by numerous influences throughout history. A genetic marker particularly suitable to investigate the paternal demographic history of populations is the non-recombining male-specific region of the Y chromosome (MSY). Using a recently established horse MSY haplotype (HT) topology and KASP™ genotyping, we illustrate MSY HT spectra of 119 Barb and Arab-Barb males, collected from the Maghreb region and European subpopulations. All detected HTs belonged to the Crown haplogroup, and the broad MSY spectrum reflects the wide variety of influential stallions throughout the breed's history. Distinct HTs and regional disparities were characterized and a remarkable number of early introduced lineages were observed. The data indicate recent refinement with Thoroughbred and Arabian patrilines, while 57% of the dataset supports historical migrations between North Africa and the Iberian Peninsula. In the Barb horse, we detected the HT linked to Godolphin Arabian, one of the Thoroughbred founders. Hence, we shed new light on the question of the ancestry of one Thoroughbred patriline. We show the strength of the horse Y chromosome as a genealogical tool, enlighten recent paternal history of North African horses, and set the foundation for future studies on the breed and the formation of conservation breeding programs.

Y-Chromosomal Insights into Breeding History and Sire Line Genealogies of Arabian Horses.

The Y chromosome is a valuable genetic marker for studying the origin and influence of paternal lineages in populations. In this study, we conducted Y-chromosomal lineage-tracing in Arabian horses. First, we resolved a Y haplotype phylogeny based on the next generation sequencing data of 157 males from several breeds. Y-chromosomal haplotypes specific for Arabian horses were inferred by genotyping a collection of 145 males representing most Arabian sire lines that are active around the globe. These lines formed three discrete haplogroups, and the same haplogroups were detected in Arabian populations native to the Middle East. The Arabian haplotypes were clearly distinct from the ones detected in Akhal Tekes, Turkoman horses, and the progeny of two Thoroughbred foundation sires. However, a haplotype introduced into the English Thoroughbred by the stallion Byerley Turk (1680), was shared among Arabians, Turkomans, and Akhal Tekes, which opens a discussion about the historic connections between Oriental horse types. Furthermore, we genetically traced Arabian sire line breeding in the Western World over the past 200 years. This confirmed a strong selection for relatively few male lineages and uncovered incongruences to written pedigree records. Overall, we demonstrate how fine-scaled Y-analysis contributes to a better understanding of the historical development of horse breeds.

Population Genetic Analysis of the Estonian Native Horse Suggests Diverse and Distinct Genetics, Ancient Origin and Contribution from Unique Patrilines.

The Estonian Native Horse (ENH) is a medium-size pony found mainly in the western islands of Estonia and is well-adapted to the harsh northern climate and poor pastures. The ancestry of the ENH is debated, including alleged claims about direct descendance from the extinct Tarpan. Here we conducted a detailed analysis of the genetic makeup and relationships of the ENH based on the genotypes of 15 autosomal short tandem repeats (STRs), 18 Y chromosomal single nucleotide polymorphisms (SNPs), mitochondrial D-loop sequence and lateral gait allele in DMRT3. The study encompassed 2890 horses of 61 breeds, including 33 ENHs. We show that the expected and observed genetic diversities of the ENH are among the highest within 52 global breeds, and the highest among 8 related Northern European ponies. The genetically closest breeds to the ENH are the Finn Horse, and the geographically more distant primitive Hucul and Konik. ENH matrilines are diverse and relate to draught and Pontic-Caspian breeds. ENH patrilines relate to draught breeds, and to a unique haplogroup not described before. None of the 33 ENHs carried the "gait" mutation, but the mutation was found in 2 Huculs. The study demonstrates that the ENH is a genetically distinct and diverse breed of ancient origin with no notable pressure of selective breeding.

The horse Y chromosome as an informative marker for tracing sire lines.

Analysis of the Y chromosome is the best-established way to reconstruct paternal family history in humans. Here, we applied fine-scaled Y-chromosomal haplotyping in horses with biallelic markers and demonstrate the potential of our approach to address the ancestry of sire lines. We de novo assembled a draft reference of the male-specific region of the Y chromosome from Illumina short reads and then screened 5.8 million basepairs for variants in 130 specimens from intensively selected and rural breeds and nine Przewalski's horses. Among domestic horses we confirmed the predominance of a young'crown haplogroup' in Central European and North American breeds. Within the crown, we distinguished 58 haplotypes based on 211 variants, forming three major haplogroups. In addition to two previously characterised haplogroups, one observed in Arabian/Coldblooded and the other in Turkoman/Thoroughbred horses, we uncovered a third haplogroup containing Iberian lines and a North African Barb Horse. In a genealogical showcase, we distinguished the patrilines of the three English Thoroughbred founder stallions and resolved a historic controversy over the parentage of the horse 'Galopin', born in 1872. We observed two nearly instantaneous radiations in the history of Central and Northern European Y-chromosomal lineages that both occurred after domestication 5,500 years ago.

Y Chromosome Uncovers the Recent Oriental Origin of Modern Stallions.

The Y chromosome directly reflects male genealogies, but the extremely low Y chromosome sequence diversity in horses has prevented the reconstruction of stallion genealogies [1, 2]. Here, we resolve the first Y chromosome genealogy of modern horses by screening 1.46 Mb of the male-specific region of the Y chromosome (MSY) in 52 horses from 21 breeds. Based on highly accurate pedigree data, we estimated the de novo mutation rate of the horse MSY and showed that various modern horse Y chromosome lineages split much later than the domestication of the species. Apart from few private northern European haplotypes, all modern horse breeds clustered together in a roughly 700-year-old haplogroup that was transmitted to Europe by the import of Oriental stallions. The Oriental horse group consisted of two major subclades: the Original Arabian lineage and the Turkoman horse lineage. We show that the English Thoroughbred MSY was derived from the Turkoman lineage and that English Thoroughbred sires are largely responsible for the predominance of this haplotype in modern horses.

High-throughput mRNA and miRNA profiling of epithelial-mesenchymal transition in MDCK cells.

BACKGROUND: Epithelial-mesenchymal transition (EMT) is an important process in embryonic development, especially during gastrulation and organ formation. Furthermore EMT is widely observed in pathological conditions, e.g., fibrosis, tumor progression and metastasis. Madin-Darby Canine Kidney (MDCK) cells are widely used for studies of EMT and epithelial plasticity. MDCK cells show an epithelial phenotype, while oncogenic Ras-transformed MDCK (MDCK-Ras) cells undergo EMT and show a mesenchymal phenotype. METHODS: RNA-Seq and miRNA-Seq analyses were performed on MDCK and MDCK-Ras cells. Data were validated by qRT-PCR. Gene signature analyses were carried out to identify pathways and gene ontology terms. For selected miRNAs target prediction was performed. RESULTS: With RNA-Seq, mRNAs of approximately half of the genes known for dog were detected. These were screened for differential regulation during Ras-induced EMT. We went further and performed gene signature analyses and found Gene Ontology (GO) terms and pathways important for epithelial polarity and implicated in EMT. Among the identified pathways, TGFß1 emerged as a central signaling factor in many EMT related pathways and biological processes. With miRNA-Seq, approximately half of the known canine miRNAs were found expressed in MDCK and MDCK-Ras cells. Furthermore, among differentially expressed miRNAs, miRNAs that are known to be important regulators of EMT were detected and new candidates were predicted. New dog miRNAs were discovered after aligning our reads to that of other species in miRBase. Importantly, we could identify 25 completely novel miRNAs with a stable hairpin structure. Two of these novel miRNAs were differentially expressed. We validated the two novel miRNAs with the highest read counts by RT-qPCR. Target prediction of a particular novel miRNA highly expressed in mesenchymal MDCK-Ras cells revealed that it targets components of epithelial cell junctional complexes. Combining target prediction for the most upregulated miRNAs and validation of the targets in MDCK-Ras cells with pathway analysis allowed us to identify two novel pathways, e.g., JAK/STAT signaling and pancreatic cancer pathways. These pathways could not be detected solely by gene set enrichment analyses of RNA-Seq data. CONCLUSION: With deep sequencing data of mRNAs and miRNAs of MDCK cells and of Ras-induced EMT in MDCK cells, differentially regulated mRNAs and miRNAs are identified. Many of the identified genes are within pathways known to be involved in EMT. Novel differentially upregulated genes in MDCK cells are interferon stimulated genes and genes involved in Slit and Netrin signaling. New pathways not yet linked to these processes were identified. A central pathway in Ras induced EMT is TGFß signaling, which leads to differential regulation of many target genes, including miRNAs. With miRNA-Seq we identified miRNAs involved in either epithelial cell biology or EMT. Finally, we describe completely novel miRNAs and their target genes.

STAT1ß is not dominant negative and is capable of contributing to gamma interferon-dependent innate immunity.

The transcription factor STAT1 is essential for interferon (IFN)-mediated immunity in humans and mice. STAT1 function is tightly regulated, and both loss- and gain-of-function mutations result in severe immune diseases. The two alternatively spliced isoforms, STAT1α and STAT1ß, differ with regard to a C-terminal transactivation domain, which is absent in STAT1ß. STAT1ß is considered to be transcriptionally inactive and to be a competitive inhibitor of STAT1α. To investigate the functions of the STAT1 isoforms in vivo, we generated mice deficient for either STAT1α or STAT1ß. As expected, the functions of STAT1α and STAT1ß in IFN-α/ß- and IFN-λ-dependent antiviral activity are largely redundant. In contrast to the current dogma, however, we found that STAT1ß is transcriptionally active in response to IFN-γ. In the absence of STAT1α, STAT1ß shows more prolonged IFN-γ-induced phosphorylation and promoter binding. Both isoforms mediate protective, IFN-γ-dependent immunity against the bacterium Listeria monocytogenes, although with remarkably different efficiencies. Our data shed new light on the potential contributions of the individual STAT1 isoforms to STAT1-dependent immune responses. Knowledge of STAT1ß's function will help fine-tune diagnostic approaches and help design more specific strategies to interfere with STAT1 activity.

Tyrosine kinase 2 promotes sepsis-associated lethality by facilitating production of interleukin-27.

The aim of this study was to test the hypothesis that gene expression and release of IL-27 may be modulated by Tyk2. Macrophages derived from the peritoneum or bone marrow of C57BL/10SnJ (WT) mice produced abundant amounts of IL-27(p28) following TLR4 activation by LPS. In contrast, production of IL-27(p28), but not EBI3, was reduced by ∼50% in TLR4-activated macrophages derived from mice with genetic deficiency of Tyk2 compared with WT macrophages. Frequencies of IL-27(p28)+F4/80+CD11b+ cells were lower in TLR4-activated macrophages derived from Tyk2-/- mice. Mechanistically, Tyk2-/- resulted in disruption of a type I IFN-dependent mechanism for production of IL-27(p28), which was induced by type I IFNs, and release of IL-27 was defective in macrophages from IFN-ß-/- and IFNAR1-/- mice. In contrast, Tyk2 was not required to mediate the effects of IL-27 on target gene expression in CD4(+) T cells. In vivo, we observed that Tyk2-/- mice have improved survival following endotoxic shock or polymicrobial sepsis induced by CLP. Plasma levels of IL-27(p28) during endotoxic shock or polymicrobial sepsis were markedly reduced in Tyk2-/- mice compared with WT mice. Disruption of IL-27 signaling using IL-27RA-/- mice was protective against sepsis-associated mortality. These data suggest that Tyk2 may mediate adverse outcomes of SIRS by promoting the production of IL-27. In conclusion, this report identifies Tyk2 as a prerequisite factor in the molecular networks that are involved in generation of IL-27.

CDK8-mediated STAT1-S727 phosphorylation restrains NK cell cytotoxicity and tumor surveillance.

The transcription factor STAT1 is important in natural killer (NK) cells, which provide immediate defense against tumor and virally infected cells. We show that mutation of a single phosphorylation site (Stat1-S727A) enhances NK cell cytotoxicity against a range of tumor cells, accompanied by increased expression of perforin and granzyme B. Stat1-S727A mice display significantly delayed disease onset in NK cell-surveilled tumor models including melanoma, leukemia, and metastasizing breast cancer. Constitutive phosphorylation of S727 depends on cyclin-dependent kinase 8 (CDK8). Inhibition of CDK8-mediated STAT1-S727 phosphorylation may thus represent a therapeutic strategy for stimulating NK cell-mediated tumor surveillance.

Organ-specific and differential requirement of TYK2 and IFNAR1 for LPS-induced iNOS expression in vivo.

Lipopolysaccharide (LPS) is an integral structural component of the outer membrane of Gram-negative bacteria and the principal active agent in the pathogenesis of endotoxin shock. LPS is a potent inducer of a variety of cytokines and inflammatory agents that lead to a profound alteration of gene expression patterns in cells and organs. The gene coding for the inducible nitric oxide synthase (iNOS) is highly responsive to LPS in vitro and in vivo and accounts for the production of nitric oxide (NO). The Janus kinase (JAK) family member tyrosine kinase 2 (TYK2) is a constituent of the interferon (IFN) type I response pathway and an important effector in the progression of endotoxin shock. Macrophages deficient for IFNalphabeta receptor chain 1 (IFNAR1) or TYK2 were shown to have an impaired LPS-induced iNOS expression. Here we determined the contribution of IFNAR1 and TYK2 to iNOS expression in vivo in a lethal LPS challenge model. TYK2 and IFNAR1 were found to be crucial for the LPS-induced iNOS mRNA and protein expression in spleen and lung that could be attributed to the Mac3-positive population. In liver LPS-induced iNOS mRNA expression was only partially impaired in TYK2-deficient mice and was unimpaired in IFNAR1-deficient mice, indicating organ specificity. TYK2(-/-) and IFNAR1(-/-) mice also differ with respect to IFNgamma production upon LPS challenge in that TYK2(-/-) mice show a defect while IFNAR1(-/-) mice do not. Our data suggest that iNOS is induced through IFNAR1 and TYK2 in Mac3-positive cells which are the main source of iNOS in spleen and lung. The LPS-induced iNOS expression in liver is independent of IFNAR1 and partially dependent on TYK2, which is most likely due to the lack of IFNgamma production in the absence of TYK2.